Near infrared adaptive optics imaging of high redshift quasars. 

Renato Falomo 

INAF - Osservatorio Astronomico di Padova, Vicolo dell'Osservatorio 5, 35122 Padova, 

Italy 

renato . f alomoSoapd . inaf . it 

Aldo Treves 

Universitd dell'Insubria, via Valleggio 11, 22100 Como, Italy 
trevesSmib . inf n . it 
Jari K. Kotilainen 

Tuorla Observatory, University of Turku, Vdisdldntie 20, FIN-21500 Piikkid, Finland 

jarkot@utu.fi 
Riccardo Scarpa 
Instituto de astrofisica de Canarias, Spain 
riccardo . scarpaOgtc . iac . es 
and 

Michela Uslenghi 
lASF-CNR Milano, Via E. Bassmi 15, Milano 1-20133, Italy 
uslenghi^mi . iasf . cnr . it 

Received ... ; accepted ... 

ABSTRACT 

The properties of high redshift quasar host galaxies are studied, in order to 
investigate the connection between galaxy evolution, nuclear activity, and the 
formation of supermassive black holes. We combine new near-IR observations 
of three high redshift quasars ( 2 < z < 3), obtained at the ESO-Very Large 
Telescope equipped with adaptive optics, with selected data from the literature. 
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For the three new objects we were able to detect and characterize the properties 
of the host galaxy, found to be consistent with those of massive elliptical galaxies 
of y\.R ~ -24.7 for the one radio loud quasar, and M^j ~ -23.8 for the two radio 
quiet quasars. When combined with existing data at lower redshift, these new 
observations depict a scenario where the host galaxies of radio loud quasars are 
seen to follow the expected trend of luminous (~5L*) elliptical galaxies under- 
going passive evolution. This trend is remarkably similar to that followed by 
radio-galaxies at z > 1.5. Radio quiet quasars hosts also follow a similar trend 
but at a lower average luminosity (~ 0.5 mag dimmer). The data indicate that 
quasar host galaxies are already fully formed at epochs as early as ~2 Gyr after 
the Big Bang and then passively fade in luminosity to the present epoch. 

Subject headings: Galaxies:active - Infrared: galaxies - Quasars: general - galaxies: 
evolution 



1. Introduction 



At low redshift quasars are hosted in o t herwise normal luminous and massive galax- 
ies jBahcall et al. 19971 : [Hamilton et al 20021 : IPunlop et al. 20031 : IPagani et al. 2003h char- 
acterized by a conspicuous spheroidal component that becomes dominant in radio loud 
objects. These galaxies appear to follow the same relationship between bulge luminosity 
and mass of th e central black hole (BH) observed in nearby inactive elliptical galaxies 
( iFerrarese 20061 . for a recent review ). If this link holds also at higher redshift the ob- 
served population of high z quasars traces the existence of ~ 10^ Mr^^ super massive BHs and 



massive spheroids at very early (< 1 Gyr) cosmic epochs (jPan et al. 200ll : iFan et al. 2003 



Willott McLure fc Jarvis 20031 ). Thi s picture seems also supported by the discovery of molec- 
ular gas and metals in high z quasars (IBertoldi et al. 20031 : iFreudling Corbin fc Korista 20031 ) , 
that are suggestive of galaxies with strong star formation. In this context it is therefore im- 
portant to push as far as possible in redshift the direct detection and characterization of 
QSO host galaxies. In particular, a key point is to probe the QSO host properties at epochs 
close to (and possibly beyond) the peak of quasar activity (z ~ 2.5). 

Until few years ago, due to the severe observational difficulties, the properties of quasar 
host galaxies at high redshift were v e ry poorly known (e.g. see the pioneering papers by 
( iHutchings 19951 : iLehnert et al. 19921 : iLowenthal et al. 19951 ). and uncertain or ambiguous 
results were produced because of inadequate quality of the images (modest resolution; low 
signal-to-noise data; non optimal analysis). 
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Deep images with adequate spatial resolution are essential. This goal is not easy to 
attain with HST because of its modest aperture that translates into a limited capabil- 
i ty to detect faint extended nebulosity unless gravitationally lensed host galaxies are used 
( iPeng et al. 20061 ). One has thus to resort to 10 meter class telescopes equipped with adap- 
tive optics (AO) systems. This keeps the advantage of both high spatial resolution and high 
sensitivity although some complications are introduced: the need for a reference point source 
close to the target, and a time and position dependent point spread function (PSF). More- 
over, unless artificial (laser) guide stars are available (not yet fully implemented in current 
AO systems) only targets which are sufficiently angularly close to relatively bright stars can 
actually be observed. 

Though the first generation of AO systems at 4m class telescopes improved the detection 
of structures of the host at low z they did not allow much improvements for distant quasars 
fautchings et al. 19981 : iHutchings et al. 19991 : iMarquez et al. 200ll : Ihacv et al. 20021 
Only very recently have AO imaging systems become available in large telescopes and could 
be used to image distant QSO with the full capability of spatial resolution and adequate 
deepness. Groom et al (2004) presented a study of 9 high z quasars imaged with the AO 
Gemini North telescope, but they were able to resolve only one radio quiet source at z= 1.93. 

In order to investigate the properties of quasar hosts at z > 2 and explore the region 
near the peak of QSO activity we are carrying out a program to secure Ks band images of 
quasars in the redshift range 2 < z < 3 using the AO system at ESO VLT. In a pre vious pilot 
work we presented the results for one radio-loud QSO (RLQ) (IFalomo et al. 20051 ). Here we 
present new observations for three high z quasars, one RLQ and two radio quiet quasars 
(RQQ). Throughout this work we use Hq = 70 km s~^ Mpc fi^ = 0.3, and = 0.7. 



Kuhlbdrodt et al 201 



2. Object selection 



Only targets sufficiently close to bright stars can be observed with adaptive optics 
systems employing natural g uide stars as reference. Because of that, we searched the latest 
( IVeron-Gettv &: Veron 20061) AGN cata l og (including data from the Sloan Digital Sky Survey 
and 2dF surveys (ISchneider et al. 20031 : iGroom et al 200ll )). for quasars in the redshift range 
2< z < 3 and 5 < 0, having a star brighter than V=14 within 30 arcsec. Beside this bright 
source needed to close the AO loop, other stars in the field of view (FoV) are necessary 
in order to characterize the PSF, both in time and position on the field of view. Thus we 
required targets to have one or more stars in the FoV for PSF characterization. Under these 
conditions the AO system at the VLT is expected to deliver images of Strehl ratio better 
than ~ 0.2 when the external seeing is <0.6 arcsec. 
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The twenty candidates fulfilling these requirements were then inspected individually 
looking at the Digitized Sky Surveys red plates. A priority was assigned according to the 
magnitude of the guide star and its distance from the target. Based on the allocated observing 
time we then chose one radio loud and 2 radio quiet objects (see Table 1). 



3. Observations and data analysis 



We acquired Ks-band images using NAOS-CONICA fiRousset et al. 20031 : iLenzen et al. 20031 ) 
the AO system on the VLT at the European Southern Observatory (ESO) in Paranal (Chile). 
The CONICA used detector was Aladdin InSb (1024x1024 pixels) that provides a field of 
view of 56x56 arcsec with a sampling of 54 mas/pixel. 

Each object was observed at random dithered positions, with small shifts applied be- 
tween successive frames, within a jitter box of ~20 arcsec around the central position of the 
object, using individual exposures of 2 minutes per frame, for a total integration time of 
38 min per observing block, each object had two observing blocks. The images (detailed in 
Table 1) were secured in service mode by ESO staff under photometric conditions. The ac- 
curacy of the photometric calibration, using standard stars observed during the same night, 
is of ±0.1 mag. 

Data reduction was performed by our own improved version of the ESO pipeline for 



jitter imaging data (IDevillard 20011 ). It first corrects for bad pixels by interpolation from 
neighboring "good" pixels, and then applies flat fielding to each image, using a normalized 
flat field obtained by subtracting and averaging a number of ON and OFF images of the 
illuminated dome. The sky background level was evaluated and sky subtraction was obtained 
for each image using an appropriate scaling and median averaging of the temporally closest 
frames. The large number of raw frames and the size of the jitter width proved to be a 
robust procedure for generating reliable sky images from the science frames themselves. All 



Table 1. Journal of the observations 



Quasar 


Type 


z 


Date 


V 


Seeing 


FWHM 


Ks 


GS(V,d) 








dd / mm / yy 


mag 


II 


II 


mag 


mag, // 


QSO 0020-304 


RQQ 


2.059 


9/12/04 


21.9 


0.5 


0.17 


17.7 


14.0, 18.5 


WGA 0633.1-233 


RLQ 


2.928 


4/02/05 


21.5 


0.6 


0.14 


18.6 


14.0, 22.3 


PKS 1041-0034 


RQQ 


2.494 


30/01/05 


20.7 


0.6 


0.15 


18.4 


14.5, 19.4 
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sky-subtracted images were then aligned to sub-pixel accuracy using 2-d cross-correlation of 
individual images using as reference all the point-like objects in the frames. 

Data for each observing block were treated separately, thus we ended up with two 
combined images for each target. These were carefully compared and found to be very 
similar, thus we further co-added them to form a single image used for all the subsequent 
analysis. 

Final modelin g of the images was done with AIDA (Astronomic al Image Decomposi- 



tion and Analysis (jUslenghi fc Falomo 20071 : lUslenghi fc Falomo 20071 )). a software package 



specifically designed to perform two dimensional model fitting of QSO images, providing 
simultaneous decomposition into the nuclear and host galaxy components. 

The most critical part of the analysis is the determination of the PSF model and the 
choice of the background level that affects the faintest external signal from the object. 



3.1. PSF modeling 

To model the PSF shape as a function of position in the FoV we used the image of the 
QSO 0633-23, that contains the largest number of stars. We first select those sources usable 
for PSF analysis on the basis of their full width at half maximum (FWHM), sharpness, 
roundness and signal-to-noise ratio, also including bright, slightly saturated stars useful to 
model the PSF faint wings. In total 14 stars were selected. 

Each star was then modeled with a gaussian for the core and an exponential function for 
the wings. Regions contaminated by close companions, saturated pixels and other evident 
defects were masked out. 

We found that both the FWHM and the ellipticity of the core component depend on 
the distance of the source from the guide star (see Figure [T]). The FWHM ranges from 0.15 
arcsec to ~ 0.3 arcsec while the ellipticity goes from 0.05 to 0.30 for objects close to the AO 
star and sources at ~40 arcsec, respectively. This analysis shows that the star most suited 
for PSF characterization should be at the same distance from the AO guide star as is the 
target. At distances from AO star <30 arcsec the size of the PSF core is stable within 10 
percent. 

We found the major axis of the PSF core is oriented along the direction connecting the 
object with the star used for the AO correction (Fig. [1]). On the contrary the shape of the 
wings is almost indepe ndent of the position in the field. This is expe cted in images obtained 
with AO corrections (( Tristram and Prieto 2005 : Cresci et al. 20051 )). 
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Fig. 1. — The first two panels show the variation of the PSF core properties with respect 
to the distance from the guide star. One pixel corresponds to 0.054 arcsec. The third 
(rightmost) panel shows the relation between direction of the elongation of the PSF core 
and the position angle of the line connecting the star with the guide star. In all panels the 
position of the target is indicated with QSO. 



A map of the PSF differences with respect to the sharpest PSF was then created by 
fitting low order polynomial curves to the PSF parameters (Figured]). This allowed us to 
evaluate the additional correction to be applied to the PSF model at the position of the 
target. The adopted PSF for the quasar was thus constructed using all the available stars 
in each frame and giving higher weight to those at similar distance from the AO star as the 
target. 

The corrections map derived using the data for QSO 0633-23 was also used for the other 
two QSOs because of the limited number of stars available in these fields. Although some 
second order variations of the PSF cannot be excluded we are confident that the general 
trends are similar. Moreover the angular separation between AO guide star and PSF star 
is similar to the one between AO star and target (difference < 5 arcsec) thus the additional 
correction is very small and cannot affect the global result of our analysis. 



3.2. QSO image decomposition 

For each source, a mask was built to exclude contamination from other sources close to 
the target, bad pixels and other possible defects. In order to take into account possible small 
variations of the local background with respect to the overall zero level of the sky subtracted 
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images we computed the average signal in a circular annulus centered on the source and 
with radii of about 3-4 and 5-6 arcsec. We checked that in such annular regions the average 
radial brightness profile of the object remains flat and that the level of the signal inside the 
annular region was consistent in all cases with the level of the background measured outside 
this annular region. This ensures that in this region there is no extra extended emission 
due to the host galaxy or associated gas. The applied correction for the level of the local 
background allows us to properly evaluate the signal of the host galaxy in the very faint 
external regions. The amount of this correction is such that only the signal below surface 
brightness fi ~ 22.5 (mag/arcsec^) is affected thus in all cases the objects would be resolved 
even without the correction. 

The QSO images were first fitted using only the point source model in order to provide 
a first check of the deviation of the target from the PSF shape. If the residuals revealed 
a significant and systematic excess over the PSF shape, the object was fitted using a two 
component model (host galaxy plus a point source). Otherwise, the object was considered 
unresolved. In all three cases presented here it was found the object was resolved thus the 
final fit of the image of the target was obtained assuming it is composed of a point source 
and an elliptical or disc galaxy convolved with the proper PSF. 

An estimate of the errors associated with the model parameters (magnitude of the nu- 
cleus, and magnitude and effective radius of the host) is shown in Fig ure El These uncertain- 



ties are consistent with those obtained using simulated quasars images (jUslenghi fc Falomo 20071 ) 



4. Results for individual objects 

J002031-3041 



This is a radio quiet quasar (z = 2.059) discovered in the 2dF survey (ICroom et al 20011 ). 
Our Ks-band image (see Figure [3]) shows the QSO close (6 arcsec) to a Ks = 13.3 star. The 
characterization of the PSF for this field is based on a bright (slightly saturated) star at 
about 6 arcsec from the target and by a couple of faint stars in the field. The core of the 
PSF was thus defined by the two faint stars while the external wing was constrained by the 
bright one (labeled PSF in Fig. [3]). The PSF was then adjusted according to the distortion 
map computed using the data of WGA J0633. 1-2333. 

The comparison of the QSO and the PSF data shows a significant residual emission up 
to ~ 1.5 arcsec from the center (see Fig. [6]) indicating that the object is resolved. The best 
2-d modeling of this object yields a host galaxy with Mr = - 23.3 and Re ~ 11 kpc. 
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Fig. 2. — The contour maps of the fit of the image of the objects as a function of 
the magnitude of the nucleus, host galaxy and the scale-length of the galaxy assuming an 
elliptical model. The three levels represent the confidence level probability of 68%, 95% and 
99% from the inner to the outer region. The dashed lines show the best fit. 



Fig. 3.— Ks-band image of the field of the QSO 0020-3041 obtained with VLT+NACO. The 
FoV shown is 40 arcsec; North is up and East to the left. The target (QSO), the guide star 
(GS) and the stars used for the characterization of the PSF are marked. 



WGA J0633. 1-2333 



This radio loud quasar (B = 21.5 ) was discovered corre lating the ROSAT WGACAT 
database with several radio catalogs ( (jPerlmann et al. 19981 )). Its redshift z = 2.928 was 
derived from prominent emission lines of Ly^+OIV, Ly a, and CIV 1540 A(Perlmann et al 
1998). The luminosity of the quasar is = -24.6. The 2-d decomposition (see Figure [6] 
) shows the host galaxy has a disturbed morphology with an extended emission structure 
at ~ 0.5 arcsec East. Assuming an elliptical model the host galaxy properties are: = 
-24.65 and Re ~ 6 kpc. 



J104117-0034 



This radio quiet quasar (V = 20.7) was discovered in the 2dF survey (jCroom et al 200ll ). 
The optical spectrum shows prominent Ly^ and CIV emission lines at z = 2.494. No radio 
emission at 20 cm is detected in the FIRST survey at the position of the quasar. 

The PSF was characterized using a star at a distance from the AO guide star similar 
to the distance between the AO guide star and the target (Fig. |5]), also adjusted using the 
map derived from the WGA J0633. 1-2333 data. Our analysis shows the object is resolved 
although we are not able to asses its morphology. Either an elliptical or a disk model for the 



Fig. 4.— NACO Ks-band image of the radio loud quasar WGA J0633.1-2333. North is up 
and east to the left. The FoV is 40 arcsec. The target (QSO), the guide star (GS) as well 
as some stars used for characterization of the PSF are marked. 

host galaxy can equally well fit the data. Under the assumption of an elliptical model the 
absolute magnitude of the host galaxy is Mji = -24.1 and the effective radius Rc ~ 2.6 kpc. 
If a disk model were assumed the host magnitude would be practically unchanged while the 
effective radius would be somewhat smaller. 

Table 2. Results of image analysis 



Quasar z Ks(nuc) Ks(host) re K-cor Re Mfl(nuc) Mfl(host) 

mag mag // Ks— >R kpc mag mag 



QSO 0020-304 


2.059 


0.33 


17.9 


19.8 


1.4 


3.05 


11.3 


-25.2 


-23.2 


QSO 0633-233 


2.928 


0.47 


19.0 


20.0 


0.7 


2.75 


5.7 


-25.0 


-24.3 


QSO 1041-003 


2.494 


0.90 


18.6 


19.6 


0.3 


2.93 


2.6 


-25.0 


-24.0 
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Fig. 5. — NACO Ks-band image of the radio quiet QSO J104117- 0034. North is up and 
east to the left. The FoV is 40 arcsec. The target (QSO), the guide star (GS) and stars used 
for characterization of the PSF are marked. 

5. Discussion 

To investigate the properties of the QSO hosts at different rcdhsifts it is preferable to 
compare data probing the same rest frame wavelengths. The Ks band at 2 < z < 3 closely 
matches the rest frame R band, thus in the following discussion our own data as well as 
data from the literature were trasformed to the R band. This was done using the cross-band 
k-correction given in Tab. 2 with details given in the Appendix. This transformation is 
virtually independent of the assumed spectral energy distribution of the host galaxy over 
the whole redshift range of interest (Am < 0.2 mag), allowing a reliable measurement of the 
rest frame luminosity. 

In order to compare our results with those published in the literature at high redshift 
we considered only observations obtained in the NIR at large (8-lOm class) telescopes or 
HST data. This allow us to perform a comparison of similar stellar populations detected in 
QSO host at lower z observed in the optical. 
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Fig. 6. — Contour plot of the object (top panel), object after subtraction of the PSF model 
(second from top), and residuals of the best fit (third from top). Scale of contour plots is 
in arcsec. Bottom panel: The radial surface brightness profile of the object (filled points) 
and the fitted model (solid line) with components (point source: dotted line; host galaxy: 
dashed line). The uncertainty in the radial surface brightness profile of the point source 
(PSF model) is about 0.1 mag at 0.5" and 0.3 mag at 1.2". 
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5.1. RLQs host galaxy evolution 



In Figure [7] we report our new measurement for the host galaxy of one RLQ at z~3, 
together with our p reviously reported RLQ hosts at z=2.55 observed with VLT+ NACO 
( iFalomo et al. 20051 ). and selected literature data. These include all results from HST 



WFPC2 images at z <0.6, our previous survey of RLQ at z<2 (IFalomo et a. 



and measurements of 4 objects by HST + NICMOS (IKukula et al. 200ll ). Two additional 



2004l : lKotilainen et al. 2007 h . 



individual points at z>2 wer e derived from the H mag of the host galaxy of lensed QSO re- 
ported by (jPeng et al. 20061 ) trasformed into R band following the method described above. 
All together these observations depict a general trend where the host luminosity increases 
by ~ 1.5 mag from present epoch up to z~3. This is fully consistent with the expected 
luminosity evolution of a massive elliptical galaxy undergoing passive evolution. On average 
this trend correspond s to that of a galaxy of luminosity ~5L* (assuming M*(Ks)=-23.9; 
( IGardner et al. 19971 ). corresponding to M*(R)=-21.1) that is undergoing passive stellar 
evolution. The dominance of an old, evolved stellar pop ulation is also supported by spectro- 
scopic studies of low redshift quasar (INolan et al. 200ll ) . 



As mentioned above it is worth to note that this result is relatively robust with respect 
to the uncertainty for the filter transformation due to the choice of SED template for the 
galaxy since the comparison is done nearly at t he same rest frarn e band. The only point 
that could move substantially is that at z = 3.27 (jPeng et al. 20061 ) since it was observed in 
H band (see Appendix). In this case if instead of an elliptical model a Sb (or Sc) SED is 
assumed the host galaxy would be 0.4 (or 0.6) magnitudes fainter suggesting a possible drop 
of the host luminosity. With the caveat of the small statistics for RLQ hosts at high redshift 
we think that the present data do not show evidence for a drop in luminosity of RLQ hosts. 



Several studies carried out at low redshift (Smith fc Heckman 1989 



Bettoni et al. 2001 



Ledlow k Owen 19951) and high redsh ift (IWillott McLure fc Jarvis 20031 : llnskip et al. 2005 



Pentericci et al 200 ll : IZirm et al 20071 ) have shown that powerful radio emission is almost 
ubiquitously linked with massive and luminous ellipticals. Indeed the global photometric and 
structural properties of radiogalaxies are identical to those of non radio early type galaxies 
of similar mass (or luminosity). This is clearly apparent at low redshift since radiog alaxies 
follow the same fundamental plane of inactive normal ellipticals (iBettoni et al. 200ll ). 



Given the above premises it is therefore of interest to compare the cosmic evolution 
of RLQ host luminosity with that of radiogalaxies (RG). (IWillott McLure &: Jarvis 20031 ) 
present a compilation of K band magnitudes of various samples of radiogalaxies. The ob- 
served K band magnitudes were converted to absolute Mr using the same transformations 
adopted for our objects. Then we have binned the data into redshift intervals of A2;=0.3 
from z=0 to z=2 and of A2;=0.5 at z > 2. The trend in luminosity of this dataset of radio- 
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galaxies is very similar to that exhibited by the host of RLQ (see Figure [7]). At z < 1 there 
is a small systematic difference (by ~ 0.5 mag) in the luminosity evolution between RLQ 
hosts and RG. This is difficult to interpret because of the r ion homogeneous definition of t he 



RG dataset (compilation from various different surveys, (IWillott McLure fc Jarvis 20031 ) ) 
can introduce selection effects in the RG samples. 

Nevertheless it is remarkable that both RG and RLQ hosts follow a similar trend of 
the luminosity up to redshift z~3. In our opinion this is suggestive of a common origin of 
the parent galaxies and also emphasizes that both types of radio loud galaxies follow the 
same evolutionary trend of inactive massive spheroids. This result can be seen also in the 
conventional K-z plot compari ng radio galaxie s and QSO hosts at z> 1 (see Figure [8]). A 



similar scenario was found by (iHutchings 20061 ) for a small sample of higher z quasars. 



5.2. RQQs host galaxy evolution 



Since the hosts of RQQ are on average less luminous than those of RLQ their study 
at high redshift is more difficult. Indeed very little is known at z>2. In Figure M we 
report our new detections for t wo radio quiet quasar hosts at z>2 compared with data frorn 
the literature at lower redshift (jDunlop et al. 20031: iKukula et al. 200ll : iHyvonen et al 20071 : 
Falomo et al. 20041 : iKotilainen et al. 2007 I : ICroom et al. 2004j l 



With the possible exception of the RQQ at z~1.9 detected by (ICroom et al. 20041 ) the 
host galaxy luminosity of RQQs appears to increase by about 1 mag from z=0 to z~ 2.5. This 
is consistent with the trend for non-active massive elliptical galaxy of M= M*-l undergoing 
simple passive evolution. This is also similar to the behavior observed for RLQs hosts but 
occur on average at a lower level of luminosity (about 0.5 mag fainter). Moreover it is worth 
to note tha t the same level of lurn inosity correspo nds to that of non- active early type galaxies 
at z ~ 1.4 Jhonghetti et al 2007h and at z ~ 1.8 flPaddi et al 2005h . 



(jPeng et al. 20061 ) report the detection of the host galaxy of a number of gravitational 
lensed quasars without radio classification, 8 in the redshift range 2 <z< 3 and 1 at z>3. 
Considering that the vast majority of QSO are radio quiet, one might assume that most of 
them are RQQ. Under this assumption and after converting their H band host luminosity to 
R band rest frame according to our SED galaxy template, we found their luminosity (<Mj^; > 
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Fig. 7. — The evolution of radio loud quasar host luminosity compared with that expected 
for massive ellipticals (at M*, M*-l and M*-2: solid, d otted and dashed line ) undergoing 
passive stellar evolution (jBressan Granato &: Silva 19981 ). The host ga laxv of the RLQ at 
z ~ 2.9 presented in this work, and another RLQ at z ~ 2.5 from (iFalomo et al. 20051 ) 
are shown as fille d pe ntagons. Other symbols re present: HST ob servations by 
( iDunlop et al. 20031) and (jPagani et al. 20031 ) (triangles), (IHooper et al. 19971) (inverted tri - 
angles); ( IKukula et al. 20011 1 (sq uares): ESQ NTT ob servat ions fjKotilainen fc Falomo 20001 ) 
(diamond). VLT observations (jFalomo et al. 20041 ) and (IKotilainen et al. 2007 I ) (circle); 
HST data for lensed hosts by ( Peng et al. 20061 ) (open stars). Each point is plotted at the 
mean redshift of the sample with error bars representing the la dispersion of the sample. In 
the case of individual objects the uncertainty of the measurement is given. A binned version 
of the data for radio galaxies shown in Fig. [8] is also given (asterisks). 
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Fig. 8. — Apparent K band magnitude versus redsh ift (the K-z plot) for luminous radio- 
galaxies (asterisks) ( (IWillott McLure &: Jarvis 20031) co mpare d with host of RLQ at z > 
1. Symbols represe nt VLT data by (jFalomo et al. 20041 ) and (IKotilainen et al. 2007 I ) (cir- 
cles): HST data b y (IKukula et al. 200ll ) (squares): HST data of grayitational lensed QSO by 
(jPeng et al. 20061 ) (open stars). The two radio loud QSO presented in this work are also 
shown (pentagons). 
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= - 24.8 at < z > ~ 2.5; M/j = -25.4 at z ~ 3.4) is well above the overall trend drawn by the 
objects at lower redshift. The only point that could be significantly affected by the choice 
of the SED is the individual source at z ^ 3.4 that, in t he case of an Sb template, would 



be fainter by 0.7 mag. In any case the (jPeng et al. 20061 ) data for these high z QSO hosts 



appear well separated from the objects at lower redshift. 

Note that the average value of Peng et al. data in the range 1 < z < 2 are fully consistent 



see Figure Ml with ou r extensive previous study performed at VLT (IFalomo et al. 2004 



Kotilainen et al. 2007 I ) thus no systematic effects due to different analysis are expected. 
One possibility to explain this difference of host average luminosity at z > 2 is that some 
(or most) of these objects are radio loud and thus their host would be naturally brighter. 
However, this may be difficult to asses because they are lensed objects. Another possibility 
is that since the objects at z > 2, studied by Peng et al., have quite luminous nuclei (average 
< Mb > = -26) also their host galaxies could be systematically more luminous than those 
in the rest of QSO samples (M^ ~ -24.3 for our RQQ at z > 2). 

Detection of extended emiss ion in the near-IR has been also found for QSO at very 



high redshift by (IHutchings 20031 ). This paper reports Gemini North direct images secured 
under average conditions (seeing 0.7 arcsec) for 5 QSO at z ~ 4.7 The observed K band host 
magnitude for these objects (as estimated from PSF removal and extrapolated fiux profile) 
are in the range Ks = 19.3 to 20.5. This corresponds to R band absolute magnitudes in 
the range M(R) = -26.6 to -27.8 assuming an elliptical galaxy SED or ~ 0.5 mag fainter 
in the case of disk galaxy SED. In both cases these data compared with those presented in 
our paper show these host galaxies are substantially more luminous (by a factor at least 5 ) 
than the trend derived from the whole dataset up to z ~ 3. However it is worth to note that 
at this redshift the K band is mapping the host galaxy at ~ 3700 A. Therefore the stellar 
population responsible for the observed emission is very different from that considered in the 
rest of the dataset (rest frame R or I). If such high luminosity will be confirmed at high z it 
may indicates a substantial amount of star formation at these epochs. 



Finally we comment on the very recent and intriguing result by ([Schramm et al 20071 ) 
of 3 RQQ at z= 2.6-2.9, that appear to have host galaxies of extremely high luminosity ( 
M^(host) = -25.8 to -26.8). They are at least 3 magnitudes above M* after including passive 
evolution. These results are very difficult to reconcile with the rest of quasar host detection 
at similar or lower redshift and would imply exceedingly high ongoing star formation. Given 
the scanty information on these luminous high z QSO it is not clear if these are exceptional 
cases possibly associated with very high luminous quasars. 

As a whole, with the caveat that at z > 2 there are only few measurements and large 
differences among objects, the present observations do not exhibit any signature of a drop 
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in luminosity (or mass) of the RQQ hosts up z~2.5. 



Conclusions 



The main motivation of this paper is to contribute to the measurement of the stellar 
luminosity of the host galaxy of high z quasars, which is a probe of the host galaxy mass. 
As a whole, while not excluding the possibility in some cases of ongoing episodes of star 
formation, the available data are consistent with no evolution in mass, indicating therefore 
that QSO host galaxies are already well formed at z ~ 3. Since then t hey passively fade to 
the present epoch. This is at odds with one of the main conclusion by (jPeng et al. 20061 ). 



This has important implications for theories of the structure formation in the Uni- 
verse. In particu lar hierarchical merging scena rios predicting a substantial mass reduction 
at early epochs (IKauffmann fc Haehnelt 20001 ). as well as those models predicting a late 
merging and assembly period for local massive spheroids, have difficulties in explaining the 
existence of a substantial population of massive, passive (red and dead) early-type galax- 
ies at high redshift (e.g. McCarthy et a l. 2004, Cimatti et al 2004, JPaddi et al 2005 



Papovich et al 2006 



(IGranato et al 2004 



Longhetti et al 20071 ) ). Only the most recent hierarchical models (e.g. 



De Lucia et al 20051 : ICroton et al 20061 : iBower et al 20061 ) which take 
into account the influence of the central supermassive black hole ( AGN feedback, e.g. 
through heating of gas in massive halos by AGN energetic), do in fact agree reasonably 
well with the observed stellar mass function and allow for the existence of massive early- 
type galaxies out to z ~4. 

Our results have important implications also for the stud y of the parameter F = ( 
Msi?/M sph), linking black hole and h ost galaxy masses (e.g. (IMerloni et al 20061 )). For 
a sample of ~ 30 quasars at z ~ 0.3 Jhabita et al. 2006h found F is consistent with the 
val ue for quiescent galaxies in the local Universe, confirming an earlier result reported 
by (IMcLure fc Dunlop 20011 ) for a different sample of low z quasars. On the other hand, 
a decreasing F was derived from a study of a sample of Seyfert galaxies at z ~ 0.36 
( iTreu et al 20071 ). At high z the situation is even le ss clear. In their study of gravitation- 
ally lensed quasar from the CASTELS-HST project (jPeng et al. 20061 ) claim F is consistent 
with the local value up to z= 1.7. After that, F sharply increases by a factor 4. Our re- 
sults strongly suggest the mass of the host galaxy does not significantly change with the 
cosmic time (at least up to z ~ 3). At face value the claim by Peng et al. would then 
imply the untenable scenario where Mbh is decreasing with the cosmic time. See however 
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Fig. 9. — The evolution of radio quiet quasar host luminosity compared with that expected 
for massive ellipticals (at M*, M*-l and M*-2; solid, dotted and da shed line, respectively ) 
undergoing passive stellar evolution (IBressan Granato fc Silva 19981 ). The two new RQQ at 



z ~ 2.05 and 2.5 presented here a re ma rked with filled circl es. The data for samples at lower 
redshift are: (jPunlop et al. 20031) and f|Pagani et al. 20031. triangles'): f|Hooper et al. 19971 ) 
( inverted triangle) ; ( Kukula et al. 200ll. squ ares): ( Hyvonen et al 20071. pentagons); 



( iFalon 



gle) 



o et al. 20041) ari d (iKotilainen et al. 2007' I) (circle); fICroom et al. 2004r) (open trian- 
(iPeng et al. 20061 ) (open stars) ; see also (IFalomo et al. 2004 ) for details on previous 
samples. Crosses represent the lumino sity of massive early typ e galaxies at z ~ 1.4 and z 



~ 1.8 studied by (ILonghetti et al 20071 ) and (jPaddi et al 20051 ). respectively. Each point is 
plotted at the mean redshift of the sample with error bars representing the la dispersion of 
the sample. In the case of individual objects the uncertainty of the measurement is given. 
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( iLauer et al. 20071) for the possi ble presence of s election bias affecting the samples consid- 
ered by (jPeng et al. 20061) and (jTreu et al 20071) . It is also worth to note the very high 
luminosities reported by (jPeng et al. 20061 ) for a number of alleged RQQ host galaxies at 
z > 2 implies an higher luminosity with respect to the passive evolution. If confirmed this 
would exacerbate the problem of the Mbh dependence on z implied by T. 

Because of the potential cosmological importance of the result in the context of the 
models of galaxy and SBH formation it is mandatory to resolve a sizeable number of objects 
at z ~ 3 and beyond. These observations should be done in the K band in order to minimize 
the uncertainty on the k-correction due to the assumed SED of the objects. 



7. Appendix - Cross filter k-correction 

We have tra sformed Ks magnitudes into rest frame R along the lines described in 
( Hogg et al 2OO2I). To perform this transformation we assumed the SED for an elliptical 



galaxy (jMannucci et al. 2001 ) and compared the integrated flux through the standard R 



(Cousins) filter to that in Ks band taking into account both the different zero point calibra- 
tion and the wavelength stretching effect by (1+z). 

In Figure [10] we report the conversion from the observed H and Ks magnitudes to 
rest frame R band for three different galaxy SED models (Elliptical, Sa and Sc) taken from 



( iMannucci et al. 20011 ). In the conversion between Ks to R the uncertainty associated to the 
choice of the SED of the galaxy is less than 0.2 mag for the whole redshift range considered 
here (from z = to z = 4). A similar uncertainty is found for the correction between H 
and R up to z = 2.5, while beyond this limit the different choice of SED lead to substantial 
different corrections up to 1 mag at z = 3.5. 
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Fig. 10. — The cross filter correction between H (bottom curves) and Ks (top curves) 
observed magnitudes at various redshift to the magnitude in R band at rest frame. Three 
different SED galaxy models are considered : Ell iptical (solid line); Sb (dashed line); Sc 
(dotted line ). Spectral templates are drawn from (jMannucci et al. 20011 1 
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